Potential impact of annual vaccination with reformulated COVID-19 vaccines: lessons from the U.S. COVID-19 Scenario Modeling Hub (preprint)

ImportanceCOVID-19 continues to cause significant hospitalizations and deaths in the United States. Its continued burden and the impact of annually reformulated vaccines remain unclear. ObjectiveTo project COVID-19 hospitalizations and deaths from April 2023-April 2025 under two plausible assumptions about immune escape (20% per year and 50% per year) and three possible CDC recommendations for the use of annually reformulated vaccines (no vaccine recommendation, vaccination for those aged 65+, vaccination for all eligible groups). DesignThe COVID-19 Scenario Modeling Hub solicited projections of COVID-19 hospitalization and deaths between April 15, 2023-April 15, 2025 under six scenarios representing the intersection of considered levels of immune escape and vaccination. State and national projections from eight modeling teams were ensembled to produce projections for each scenario. SettingThe entire United States. ParticipantsNone. ExposureAnnually reformulated vaccines assumed to be 65% effective against strains circulating on June 15 of each year and to become available on September 1. Age and state specific coverage in recommended groups was assumed to match that seen for the first (fall 2021) COVID-19 booster. Main outcomes and measuresEnsemble estimates of weekly and cumulative COVID-19 hospitalizations and deaths. Expected relative and absolute reductions in hospitalizations and deaths due to vaccination over the projection period. ResultsFrom April 15, 2023-April 15, 2025, COVID-19 is projected to cause annual epidemics peaking November-January. In the most pessimistic scenario (high immune escape, no vaccination recommendation), we project 2.1 million (90% PI: 1,438,000-4,270,000) hospitalizations and 209,000 (90% PI: 139,000-461,000) deaths, exceeding pre-pandemic mortality of influenza and pneumonia. In high immune escape scenarios, vaccination of those aged 65+ results in 230,000 (95% CI: 104,000-355,000) fewer hospitalizations and 33,000 (95% CI: 12,000-54,000) fewer deaths, while vaccination of all eligible individuals results in 431,000 (95% CI: 264,000-598,000) fewer hospitalizations and 49,000 (95% CI: 29,000-69,000) fewer deaths. Conclusion and RelevanceCOVID-19 is projected to be a significant public health threat over the coming two years. Broad vaccination has the potential to substantially reduce the burden of this disease. Key pointsO_ST_ABSQuestionC_ST_ABSWhat is the likely impact of COVID-19 from April 2023-April 2025 and to what extent can vaccination reduce hospitalizations and deaths? FindingsUnder plausible assumptions about viral evolution and waning immunity, COVID-19 will likely cause annual epidemics peaking in November-January over the two-year projection period. Though significant, hospitalizations and deaths are unlikely to reach levels seen in previous winters. The projected health impacts of COVID-19 are reduced by 10-20% through moderate use of reformulated vaccines. MeaningCOVID-19 is projected to remain a significant public health threat. Annual vaccination can reduce morbidity, mortality, and strain on health systems.

Informing pandemic response in the face of uncertainty. An evaluation of the U.S. COVID-19 Scenario Modeling Hub (preprint)

Our ability to forecast epidemics more than a few weeks into the future is constrained by the complexity of disease systems, our limited ability to measure the current state of an epidemic, and uncertainties in how human action will affect transmission. Realistic longer-term projections (spanning more than a few weeks) may, however, be possible under defined scenarios that specify the future state of critical epidemic drivers, with the additional benefit that such scenarios can be used to anticipate the comparative effect of control measures. Since December 2020, the U.S. COVID-19 Scenario Modeling Hub (SMH) has convened multiple modeling teams to make 6-month ahead projections of the number of SARS-CoV-2 cases, hospitalizations and deaths. The SMH released nearly 1.8 million national and state-level projections between February 2021 and November 2022. SMH performance varied widely as a function of both scenario validity and model calibration. Scenario assumptions were periodically invalidated by the arrival of unanticipated SARS-CoV-2 variants, but SMH still provided projections on average 22 weeks before changes in assumptions (such as virus transmissibility) invalidated scenarios and their corresponding projections. During these periods, before emergence of a novel variant, a linear opinion pool ensemble of contributed models was consistently more reliable than any single model, and projection interval coverage was near target levels for the most plausible scenarios (e.g., 79% coverage for 95% projection interval). SMH projections were used operationally to guide planning and policy at different stages of the pandemic, illustrating the value of the hub approach for long-term scenario projections.

Challenges of COVID-19 Case Forecasting in the US, 2020-2021 (preprint)

During the COVID-19 pandemic, forecasting COVID-19 trends to support planning and response was a priority for scientists and decision makers alike. In the United States, COVID-19 forecasting was coordinated by a large group of universities, companies, and government entities led by the Centers for Disease Control and Prevention and the US COVID-19 Forecast Hub (https://covid19forecasthub.org). We evaluated approximately 9.7 million forecasts of weekly state-level COVID-19 cases for predictions 1-4 weeks into the future submitted by 24 teams from August 2020 to December 2021. We assessed coverage of central prediction intervals and weighted interval scores (WIS), adjusting for missing forecasts relative to a baseline forecast, and used a Gaussian generalized estimating equation (GEE) model to evaluate differences in skill across epidemic phases that were defined by the effective reproduction number. Overall, we found high variation in skill across individual models, with ensemble-based forecasts outperforming other approaches. Forecast skill relative to the baseline was generally higher for larger jurisdictions (e.g., states compared to counties). Over time, forecasts generally performed worst in periods of rapid changes in reported cases (either in increasing or decreasing epidemic phases) with 95% prediction interval coverage dropping below 50% during the growth phases of the winter 2020, Delta, and Omicron waves. Ideally, case forecasts could serve as a leading indicator of changes in transmission dynamics. However, while most COVID-19 case forecasts outperformed a naive baseline model, even the most accurate case forecasts were unreliable in key phases. Further research could improve forecasts of leading indicators, like COVID-19 cases, by leveraging additional real-time data, addressing performance across phases, improving the characterization of forecast confidence, and ensuring that forecasts were coherent across spatial scales. In the meantime, it is critical for forecast users to appreciate current limitations and use a broad set of indicators to inform pandemic-related decision making. Author SummaryAs SARS-CoV-2 began to spread throughout the world in early 2020, modelers played a critical role in predicting how the epidemic could take shape. Short-term forecasts of epidemic outcomes (for example, infections, cases, hospitalizations, or deaths) provided useful information to support pandemic planning, resource allocation, and intervention. Yet, infectious disease forecasting is still a nascent science, and the reliability of different types of forecasts is unclear. We retrospectively evaluated COVID-19 case forecasts, which were often unreliable. For example, forecasts did not anticipate the speed of increase in cases in early winter 2020. This analysis provides insights on specific problems that could be addressed in future research to improve forecasts and their use. Identifying the strengths and weaknesses of forecasts is critical to improving forecasting for current and future public health responses.

Heterogeneous changes in mobility in response to the SARS-CoV-2 Omicron BA.2 outbreak in Shanghai (preprint)

The coronavirus disease 2019 (COVID-19) pandemic and the measures taken by authorities to control its spread had altered human behavior and mobility patterns in an unprecedented way. However, it remains unclear whether the population response to a COVID-19 outbreak varies within a city or among demographic groups. Here we utilized passively recorded cellular signaling data at a spatial resolution of 1km x 1km for over 5 million users and epidemiological surveillance data collected during the SARS-CoV-2 Omicron BA.2 outbreak from February to June 2022 in Shanghai, China, to investigate the heterogeneous response of different segments of the population at the within-city level and examine its relationship with the actual risk of infection. Changes in behavior were spatially heterogenous within the city and population groups, and associated with both the infection incidence and adopted interventions. We also found that males and individuals aged 30-59 years old traveled more frequently, traveled longer distances, and their communities were more connected; the same groups were also associated with the highest SARS-CoV-2 incidence. Our results highlight the heterogeneous behavioral change of the Shanghai population to the SARS-CoV-2 Omicron BA.2 outbreak and the its effect on the heterogenous spread of COVID-19, both spatially and demographically. These findings could be instrumental for the design of targeted interventions for the control and mitigation of future outbreaks of COVID-19 and, more broadly, of respiratory pathogens.

Spatial scales of COVID-19 transmission in Mexico (preprint)

During outbreaks of emerging infectious diseases, internationally connected cities often experience large and early outbreaks, while rural regions follow after some delay. This hierarchical structure of disease spread is influenced primarily by the multiscale structure of human mobility. However, during the COVID-19 epidemic, public health responses typically did not take into consideration the explicit spatial structure of human mobility when designing non-pharmaceutical interventions (NPIs). NPIs were applied primarily at national or regional scales. Here we use weekly anonymized and aggregated human mobility data and spatially highly resolved data on COVID-19 cases, deaths and hospitalizations at the municipality level in Mexico to investigate how behavioural changes in response to the pandemic have altered the spatial scales of transmission and interventions during its first wave (March - June 2020). We find that the epidemic dynamics in Mexico were initially driven by SARS-CoV-2 exports from Mexico State and Mexico City, where early outbreaks occurred. The mobility network shifted after the implementation of interventions in late March 2020, and the mobility network communities became more disjointed while epidemics in these communities became increasingly synchronised. Our results provide actionable and dynamic insights into how to use network science and epidemiological modelling to inform the spatial scale at which interventions are most impactful in mitigating the spread of COVID-19 and infectious diseases in general.

Characterizing collective physical distancing in the U.S. during the first nine months of the COVID-19 pandemic (preprint)

The COVID-19 pandemic offers an unprecedented natural experiment providing insights into the emergence of collective behavioral changes of both exogenous (government mandated) and endogenous (spontaneous reaction to infection risks) origin. Here, we characterize collective physical distancing -- mobility reductions, minimization of contacts, shortening of contact duration -- in response to the COVID-19 pandemic in the pre-vaccine era by analyzing de-identified, privacy-preserving location data for a panel of over 5.5 million anonymized, opted-in U.S. devices. We define five indicators of users' mobility and proximity to investigate how the emerging collective behavior deviates from the typical pre-pandemic patterns during the first nine months of the COVID-19 pandemic. We analyze both the dramatic changes due to the government mandated mitigation policies and the more spontaneous societal adaptation into a new (physically distanced) normal in the fall 2020. The indicators defined here allow the quantification of behavior changes across the rural/urban divide and highlight the statistical association of mobility and proximity indicators with metrics characterizing the pandemic's social and public health impact such as unemployment and deaths. This study provides a framework to study massive social distancing phenomena with potential uses in analyzing and monitoring the effects of pandemic mitigation plans at the national and international level.

The unequal effects of the health-economy tradeoff during the COVID-19 pandemic (preprint)

The potential tradeoff between health outcomes and economic impact has been a major challenge in the policy making process during the COVID-19 pandemic. Epidemic-economic models designed to address this issue are either too aggregate to consider heterogeneous outcomes across socio-economic groups, or, when sufficiently fine-grained, not well grounded by empirical data. To fill this gap, we introduce a data-driven, granular, agent-based model that simulates epidemic and economic outcomes across industries, occupations, and income levels with geographic realism. The key mechanism coupling the epidemic and economic modules is the reduction in consumption demand due to fear of infection. We calibrate the model to the first wave of COVID-19 in the New York metropolitan area, showing that it reproduces key epidemic and economic statistics, and then examine counterfactual scenarios. We find that: (a) both high fear of infection and strict restrictions similarly harm the economy but reduce infections; (b) low-income workers bear the brunt of both the economic and epidemic harm; (c) closing non-customer-facing industries such as manufacturing and construction only marginally reduces the death toll while considerably increasing unemployment; and (d) delaying the start of protective measures does little to help the economy and worsens epidemic outcomes in all scenarios. We anticipate that our model will help designing effective and equitable non-pharmaceutical interventions that minimize disruptions in the face of a novel pandemic.

Estimating the impact of COVID-19 vaccine allocation inequities: a modeling study (preprint)

Access to COVID-19 vaccines on the global scale has been drastically impacted by structural socio-economic inequities. Here, we develop a data-driven, age-stratified epidemic model to evaluate the effects of COVID-19 vaccine inequities in twenty lower middle and low income countries (LMIC) sampled from all WHO regions. We focus on the first critical months of vaccine distribution and administration, exploring counterfactual scenarios where we assume the same per capita daily vaccination rate reported in selected high income countries. We estimate that, in this high vaccine availability scenario, more than 50% of deaths (min-max range: 56%-99%) that occurred in the analyzed countries could have been averted. We further consider a scenario where LMIC had similarly early access to vaccine doses as high income countries; even without increasing the number of doses, we estimate an important fraction of deaths (min-max range: 7%-73%) could have been averted. In the absence of equitable allocation, the model suggests that considerable additional non-pharmaceutical interventions would have been required to offset the lack of vaccines (min-max range: 15%-75%). Overall, our results quantify the negative impacts of vaccines inequities and call for amplified global efforts to provide better access to vaccine programs in low and lower middle income countries

Forecasting hospital-level COVID-19 admissions using real-time mobility data (preprint)

For each of the COVID-19 pandemic waves, hospitals have had to plan for deploying surge capacity and resources to manage large but transient increases in COVID-19 admissions. While a lot of effort has gone into predicting regional trends in COVID-19 cases and hospitalizations, there are far fewer successful tools for creating accurate hospital-level forecasts. At the same time, anonymized phone-collected mobility data proved to correlate well with the number of cases for the first two waves of the pandemic (spring 2020, and fall-winter 2021). In this work, we show how mobility data could bolster hospital-specific COVID-19 admission forecasts for five hospitals in Massachusetts during the initial COVID-19 surge. The high predictive capability of the model was achieved by combining anonymized, aggregated mobile device data about users' contact patterns, commuting volume, and mobility range with COVID hospitalizations and test-positivity data. We conclude that mobility-informed forecasting models can increase the lead-time of accurate predictions for individual hospitals, giving managers valuable time to strategize how best to allocate resources to manage forthcoming surges.

Impact of SARS-CoV-2 vaccination of children ages 5-11 years on COVID-19 disease burden and resilience to new variants in the United States, November 2021-March 2022: a multi-model study (preprint)

Background: SARS-CoV-2 vaccination of persons aged 12 years and older has reduced disease burden in the United States. The COVID-19 Scenario Modeling Hub convened multiple modeling teams in September 2021 to project the impact of expanding vaccine administration to children 5-11 years old on anticipated COVID-19 burden and resilience against variant strains. Methods: Nine modeling teams contributed state- and national-level projections for weekly counts of cases, hospitalizations, and deaths in the United States for the period September 12, 2021 to March 12, 2022. Four scenarios covered all combinations of: 1) presence vs. absence of vaccination of children ages 5-11 years starting on November 1, 2021; and 2) continued dominance of the Delta variant vs. emergence of a hypothetical more transmissible variant on November 15, 2021. Individual team projections were combined using linear pooling. The effect of childhood vaccination on overall and age-specific outcomes was estimated by meta-analysis approaches. Findings: Absent a new variant, COVID-19 cases, hospitalizations, and deaths among all ages were projected to decrease nationally through mid-March 2022. Under a set of specific assumptions, models projected that vaccination of children 5-11 years old was associated with reductions in all-age cumulative cases (7.2%, mean incidence ratio [IR] 0.928, 95% confidence interval [CI] 0.880-0.977), hospitalizations (8.7%, mean IR 0.913, 95% CI 0.834-0.992), and deaths (9.2%, mean IR 0.908, 95% CI 0.797-1.020) compared with scenarios where children were not vaccinated. This projected effect of vaccinating children 5-11 years old increased in the presence of a more transmissible variant, assuming no change in vaccine effectiveness by variant. Larger relative reductions in cumulative cases, hospitalizations, and deaths were observed for children than for the entire U.S. population. Substantial state-level variation was projected in epidemic trajectories, vaccine benefits, and variant impacts. Conclusions: Results from this multi-model aggregation study suggest that, under a specific set of scenario assumptions, expanding vaccination to children 5-11 years old would provide measurable direct benefits to this age group and indirect benefits to the all-age U.S. population, including resilience to more transmissible variants.

Preliminary modeling estimates of the relative transmissibility and immune escape of the Omicron SARS-CoV-2 variant of concern in South Africa (preprint)

We develop a stochastic, multi-strain, compartmental epidemic model to estimate the relative transmissibility and immune escape of the Omicron variant of concern (VOC) in South Africa. The model integrates population, non-pharmaceutical interventions, vaccines, and epidemiological data and it is calibrated in the period May 1st, 2021 - November 23rd, 2021. We explore a parameter space of relative transmissibility with respect to the Delta variant and immune escape for Omicron by assuming an initial seeding, from unknown origin, in the first week of October 2021. We identify a region of the parameter space where combinations of relative transmissibility and immune escape are compatible with the growth of the epidemic wave. We also find that changes in the generation time associated with Omicron infections strongly affect the results concerning its relative transmissibility. The presented results are informed by current knowledge of Omicron and subject to changes.

Anatomy of the first six months of COVID-19 Vaccination Campaign in Italy (preprint)

We analyze the effectiveness of the first six months of vaccination campaign against SARS-CoV-2 in Italy by using a computational epidemic model which takes into account demographic, mobility, vaccines, as well as estimates of the introduction and spreading of the more transmissible Alpha variant. We consider six sub-national regions and study the effect of vaccines in terms of number of averted deaths, infections, and reduction in the Infection Fatality Rate (IFR) with respect to counterfactual scenarios with the actual non-pharmaceuticals interventions but no vaccine administration. Furthermore, we compare the effectiveness in counterfactual scenarios with different vaccines allocation strategies and vaccination rates. Our results show that, as of 2021/07/05, vaccines averted 29, 350 ( IQR [16, 454 − 42, 826]) deaths and 4, 256, 332 ( IQR [1, 675, 564 − 6, 980, 070]) infections and a new pandemic wave in the country. During the same period, they achieved a −22.2% ( IQR [−31.4%; −13.9%]) reduction in the IFR. We show that a campaign that would have strictly prioritized age groups at higher risk of dying from COVID-19, besides frontline workers, would have implied additional benefits both in terms of avoided fatalities and reduction in the IFR. Strategies targeting the most active age groups would have prevented a higher number of infections but would have been associated with more deaths. Finally, we study the effects of different vaccination intake scenarios by rescaling the number of available doses in the time period under study to those administered in other countries of reference. The modeling framework can be applied to other countries to provide a mechanistic characterization of vaccination campaigns worldwide.

Higher education responses to COVID-19 in the United States: Evidence for the impacts of university policy (preprint)

With a dataset of testing and case counts from over 1,400 institutions of higher education (IHEs) in the United States, we analyze the number of infections and deaths from SARS-CoV-2 in the counties surrounding these IHEs during the Fall 2020 semester (August to December, 2020). We used a matching procedure designed to create groups of counties that are aligned along age, race, income, population, and urban/rural categories---socio-demographic variables that have been shown to be correlated with COVID-19 outcomes. We find that counties with IHEs that remained primarily online experienced fewer cases and deaths during the Fall 2020 semester; whereas before and after the semester, these two groups had almost identical COVID-19 incidence. Additionally, we see fewer deaths in counties with IHEs that reported conducting any on-campus testing compared to those that reported none. We complement the statistical analysis with a case study of IHEs in Massachusetts---a rich data state in our dataset---which further highlights the importance of IHE-affiliated testing for the broader community. The results in this work suggest that campus testing can itself be thought of as a mitigation policy and that allocating additional resources to IHEs to support efforts to regularly test students and staff would be beneficial to mitigating the spread of COVID-19 in the general population.

Projected resurgence of COVID-19 in the United States in July-December 2021 resulting from the increased transmissibility of the Delta variant and faltering vaccination (preprint)

What is already known about this topic?The highly transmissible SARS-CoV-2 Delta variant has begun to cause increases in cases, hospitalizations, and deaths in parts of the United States. With slowed vaccination uptake, this novel variant is expected to increase the risk of pandemic resurgence in the US in July--December 2021. What is added by this report?Data from nine mechanistic models project substantial resurgences of COVID-19 across the US resulting from the more transmissible Delta variant. These resurgences, which have now been observed in most states, were projected to occur across most of the US, coinciding with school and business reopening. Reaching higher vaccine coverage in July--December 2021 reduces the size and duration of the projected resurgence substantially. The expected impact of the outbreak is largely concentrated in a subset of states with lower vaccination coverage. What are the implications for public health practice?Renewed efforts to increase vaccination uptake are critical to limiting transmission and disease, particularly in states with lower current vaccination coverage. Reaching higher vaccination goals in the coming months can potentially avert 1.5 million cases and 21,000 deaths and improve the ability to safely resume social contacts, and educational and business activities. Continued or renewed non-pharmaceutical interventions, including masking, can also help limit transmission, particularly as schools and businesses reopen.

Model-based evaluation of alternative reactive class closure strategies against COVID-19 (preprint)

There are contrasting results concerning the effect of reactive school closure on SARS-CoV-2 transmission. To shed light on this controversy, here we develop a data-driven computational model of SARS-CoV-2 transmission to investigate mechanistically the effect on COVID-19 outbreaks of school closure strategies based on syndromic surveillance and antigen screening of students. We found that by reactively closing classes based on syndromic surveillance, SARS-CoV-2 infections are reduced by no more than 13.1% (95%CI: 8.6%-20.2 %), due to the low probability of timely symptomatic case identification among the young population. We thus investigated an alternative triggering mechanism based on repeated screening of students using antigen tests. Should population-level social distancing measures unrelated to schools enable maintaining the reproduction number (R) at 1.3 or lower, an antigen-based screening strategy is estimated to fully prevent COVID-19 outbreaks in the general population. Depending on the contribution of schools to transmission, this strategy can either prevent COVID-19 outbreaks for R up to 1.9 or to at least greatly reduce outbreak size in very conservative scenarios about school contribution to transmission. Moving forward, the adoption of antigen-based screenings in schools could be instrumental to limit COVID-19 burden while vaccines continue to roll out through 2021, especially in light of possible continued emergence of SARS-CoV-2 variants.

Cryptic transmission of SARS-CoV-2 and the first COVID-19 wave in Europe and the United States (preprint)

Given the narrowness of the initial testing criteria, the SARS-CoV-2 virus spread through cryptic transmission in January and February, setting the stage for the epidemic wave experienced in March and April, 2020. We use a global metapopulation epidemic model to provide a mechanistic understanding of the global dynamic underlying the establishment of the COVID-19 pandemic in Europe and the United States (US). The model is calibrated on international case introductions at the early stage of the pandemic. We find that widespread community transmission of SARS-CoV-2 was likely in several areas of Europe and the US by January 2020, and estimate that by early March, only 1 - 3 in 100 SARS-CoV-2 infections were detected by surveillance systems. Modeling results indicate international travel as the key driver of the introduction of SARS-CoV-2 with possible importation and transmission events as early as December, 2019. We characterize the resulting heterogeneous spatio-temporal spread of SARS-CoV-2 and the burden of the first COVID-19 wave (February-July 2020). We estimate infection attack rates ranging from 0.78%-15.2% in the US and 0.19%-13.2% in Europe. The spatial modeling of SARS-CoV-2 introductions and spreading provides insights into the design of innovative, model-driven surveillance systems and preparedness plans that have a broader initial capacity and indication for testing.

Data-driven estimate of SARS-CoV-2 herd immunity threshold in populations with individual contact pattern variations (preprint)

The development of efficacious vaccines has made it possible to envision mass vaccination programs aimed at suppressing SARS-CoV-2 transmission around the world. Here we use a data-driven age-structured multilayer representation of the population of 34 countries to estimate the disease induced immunity threshold, accounting for the contact variability across individuals. We show that the herd immunization threshold of random (un-prioritized) mass vaccination programs is generally larger than the disease induced immunity threshold. We use the model to test two additional vaccine prioritization strategies, transmission-focused and age-based, in which individuals are inoculated either according to their behavior (number of contacts) or infection fatality risk, respectively. Our results show that in the case of a sterilizing vaccine the behavioral strategy achieves herd-immunity at a coverage comparable to the disease-induced immunity threshold, but it appears to have inferior performance in averting deaths than the risk vaccination strategy. The presented results have potential use in defining the effects that the heterogeneity of social mixing and contact patterns has on herd immunity levels and the deployment of vaccine prioritization strategies.

Estimating the spreading and dominance of SARS-CoV-2 VOC 202012/01 (lineage B.1.1.7) across Europe (preprint)

We develop a two strain, age-structured, compartmental model to assess the spreading potential of the B.1.1.7 variant across several European metropolitan areas and countries. The model accounts for B.1.1.7 introductions from the UK and different locations, as well as local mitigation policies in the time period 2020/09 - 2021/02. In the case of an increase of transmissibility of 50%, the B.1.1.7 variant has the potential to become dominant in all investigated areas by the end of March 2021.

DeepGLEAM: A hybrid mechanistic and deep learning model for COVID-19 forecasting (preprint)

We introduce DeepGLEAM, a hybrid model for COVID-19 forecasting. DeepGLEAM combines a mechanistic stochastic simulation model GLEAM with deep learning. It uses deep learning to learn the correction terms from GLEAM, which leads to improved performance. We further integrate various uncertainty quantification methods to generate confidence intervals. We demonstrate DeepGLEAM on real-world COVID-19 mortality forecasting tasks.